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If  in  a stratified  ocean  an  acoustical  signal  is  emitted  by  a source  at  a known 
time  and  received  at  another  location  at  a known  time,  if  the  angle  of  arrival 
of  the  ray  path  at  the  receiver  is  known,  and  if  the  speed  of  sound  profile  is 
known,  then  the  acoustical  path  from  receiver  to  source  can  be  traced  back  and 
the  position  of  the  source  estimated. 

The  relative  accuracies  of  two  models  of  speed-of-sound  profiles  used  in  per- 
forming this  task  are  compared.  The  discrepancies  in  trace-back  between  a 
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20.  speed  of  sound  profile  modeled  by  a succession  of  isospeed  layers  and 

modeled  by  a succession  of  isogradient  layers  are  obtained  for  two  situ- 
ations: (1)  when  the  value  of  the  speed  of  sound  in  each  isospeed  layer 
is  the  average  over  depth  of  the  speed  of  sound  in  the  analogous  isogradien 
layer,  and  (2)  when  the  angle  of  elevation  of  the  ray  in  the  isospeed 
layer  is  the  average  of  the  entrance  and  exit  angles  of  the  equivalent  ray 
in  the  analogous  isogradient  layer.  It  is  found  that  the  second  situation 
leads  to  discrepancies  in  range  and  depth  which  are  always  smaller  than 
those  encountered  in  the  first,  and  in  particular  the  discrepancies  in 
the  first  situation  are  orders  of  magnitude  smaller  than  those  in  the  sec- 
ond for  cases  wherein  the  rays  tend  to  vertex. 

The  importance  of  the  Snell's  law  invariant  (cos  0)/c(z)  is  pointed  out. 
Choosing  an  incorrect  value  for  either  the  angle  0 or  the  actual  speed  of 
sound  c(z)  at  the  receiver  can  lead  to  errors  in  localization  of  the 
source  which  are  unacceptably  large  for  accurate  trace-back  in  situations 
where  the  ray  is  nearly  horizontal  over  some  part  of  its  path. 
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I . INTRODUCTION 


Assume  that  an  acoustical  source  initiates  a signal  at  a known  time  and 
that  this  signal  is  subsequently  received  by  a directional  array.  If  the 
time  of  arrival  of  the  signal  at  each  hydrophone  of  the  array  is  known, 
and  the  speed-of -sound  profile  of  the  water  column  is  known,  it  is  then 
possible  to  determine  the  angle  9r  of  the  arrival  of  the  signal  at  the 
array  and  trace  the  trajectory  of  the  signal  back  over  the  relevant  ray 
path.  With  the  time  of  flight  known,  the  original  position  (range  and 
depth)  of  the  acoustical  source  can  then  be  estimated. 

Of  the  various  sources  of  error  possible  in  implementing  this  process , 
this  report  is  directed  toward  those  arising  from  the  fact  that  the  speed- 
of-sound  profile  used  in  calculating  the  ray  path  is  modeled  on  the  available 
experimental  data  sampling  the  true  profile. 

If  the  water  column  is  divided  into  horizontal  layers,  the  speed-of-sound 
profile  in  each  of  these  layers  can  be  approximated  with  some  simple  function 
of  depth  which  models  the  observed  profile  and  allows  rapid  calculation  of 
the  ray  path  and  time  of  flight  through  that  layer. 

The  simplest  model  would  be  a one -parameter  fit:  the  assumption  of  an 
isospeed  layer  whose  speed  of  sound  c^  is  obtained  from  analysis  of  the  data 
within  that  layer.  The  ray  path  would  be  a straight  line  whose  angle  of 
elevation  would  be  determined  by  a layer-by-layer  application  of  Snell's  Law 
along  the  ray  path.  This  model  is  suggested  in  Figure  1.1a. 

A more  sophisticated  model  would  be  a two-parameter  fit:  the  assumption 
of  an  isogradient  layer  whose  speed  of  sound  c(z)  is  a linear  function  of 
depth  z , 

c(z)  = c (z  ) + g(z  - z ) 
o o 
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where  g is  the  gradient  and  zq  is  some  reference  depth  in  the  layer.  Here, 
the  two  parameters  c(zq)  and  g would  be  obtained  by  fitting  c(z)  to  the  data 
representing  the  observed  profile  within  the  layer.  The  ray  path  in  this 
case  would  be  a circular  arc  of  fixed  radius  within  each  layer.  See  Figure 
1.1b,  and  notice  in  particular  that  the  two  parameter  fit  allows  c to  be 
discontinuous  across  the  top  or  bottom  of  a layer. 

A commonly-encountered  isogradient-layer  model  is  the  "constrained"  two 
parameter  fit  (Fig.  1.1c),  for  which  the  speed  of  sound  is  not  allowed  to 
have  discontinuities.  A moment's  thought  reveals  that  this  constraint  must 
reduce  the  accuracy  of  this  model  below  that  of  the  (unconstrained)  two 
parameter  fit  of  Figure  1.1b. 

Other  models  with  more  parameters  can  be  postulated  (corresponding  to 
fitting  the  data  with  power  series  in  z of  increasing  order) , but  the 
resulting  ray  paths  become  more  complicated  than  usually  desired. 

Given  that  the  experimental  data  are  accurate,  models  with  more  parameters 
would  be  better  capable  of  reproducing  the  experimentally  determined  profile, 
at  least  until  there  cure  too  many  parameters  for  the  number  of  data  or  until 
whatever  desired  smoothing  of  the  data  is  lost.  Below  these  limits,  however, 
the  two-parameter  fit  should  produce  a ray  path  which  conforms  more  closely 
to  that  in  "real"  profile  than  should  the  one-parameter  fit.  It  is  plausible 
that  this  should  also  be  true  for  the  calculations  of  the  time  of  flight 
through  the  layers. 

In  this  report  we  compare  the  predictions  of  the  isogradient  model  with 
those  of  two  isospeed  models . Our  approach  will  be  to  take  the  isogradient 
model  as  the  reference  and  compare  the  predictions  of  two  isospeed  models 
with  those  of  the  isogradient  model.  That  isospeed  model  whose  predictions 
conform  more  closely  to  the  reference  will  be  deemed  the  better  of  the  two 
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with  respect  to  locating  the  acoustical  source.  We  also  consider  some 


effects  on  accurate  prediction  arising  from  an  incorrect  value  of  either 
the  angle  of  arrival  6^  or  the  speed  of  sound  c at  the  array. 


t 
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II.  FORMULATION 


Both  isogradient-layer  models  and  isospeed-layer  models  approximate 
the  true  speed-of-sound  profile  for  the  purpose  of  tracing  received 
acoustical  signals  back  to  their  apparent  source. 

What  is  known  about  a received  signal  for  the  situation  considered 
here  is  the  time-of-flight  of  the  signal  from  source  to  receiver  and  the 
angle  6 with  which  the  signal  impinges  on  the  receiver.  This  angle, 
together  with  prior  information  about  the  speed  of  sound  at  the  receiver, 
cr>  provides  a specific  value  for  the  ray  invariant  (cos  9r/cr)  which  uniquely 
labels  the  ray.  It  follows  that  comparison  of  isogradient  and  isospeed 
models  must  be  based  on  studying  the  different  ray  paths  and  travel  times 
resulting  from  tracing  back  a ray  of  specified  invariant . 

A.  Isogradient  Layer 

Refer  to  Fig.  1(a).  Let  a ray  enter  an  isogradient  layer  of  thickness 
Az  at  an  angle  of  elevation  9 and  leave  the  layer  with  an  angle  of  eleva- 
tion 9 + A9 . The  gradient  g is  defined  so  that  positive-gradient  water 
has  speed  of  sound  c(z)  increasing  with  increasing  depth,  and  the  speed  of 
sound  at  the  bottom  of  the  layer  is  c.  It  is  useful  to  define  a quantity 

D = & Az  2.1 

c 

which  is  the  fractional  change  of  the  speed  of  sound  Ac/c  across  the  layer. 
Application  of  Snell's  Law  along  the  ray  relates  D and  A9, 

cos  9 _ cos  (9  + A9)  ^ 

c = c (1  - D) 


Use  of  cos  (0  + A9)  = cos  0 cos  A0  - sin  G sin  A9,  assuming  that  the 
incremental  angle  of  A9  is  very  small  so  that  cos  A0  and  sin  A9  can  be 
replaced  by  truncated  power  series  in  A0,  and  collection  of  like-order 
terms  gives 

D = tan  0 A0  +(l/2)(A0)2  - (l/6)  tan  0 (A0)3  - (l/24)  (A0)4  + ...  2.3 

A0  « 1 . 


Between  entering  and  leaving  the  layer  the  range  has  increased  by 

A^i.  , and  the  time-of-flight  has  increased  by  At  . Since  the  ray  path  is 
g g 

a section  of  a circular  arc,  elementary  analysis  yields  the  known  expressions 
= g cos"  9 ^sin  (e  + A0)  - sin  2-4 


and 


At  = - In 
g 


1_ 

g 


sin  1.0 


c(l  - D)  1 + sin  0 


an  ) 
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In  the  first  expression,  Eq.  2.4,  expansion  of  the  first  sine  and  use  of 
A0  <<  1 results  in  the  power  series 

A K = j [A0  -(1/2)  tan  9(A0)2  - (l/6)(A9)  3 + (l/24)  tan  0 (A0)4  +...]  . 2.6 

o © 

Now,  use  of  Eqs . 2.1  and  2.3  yields  the  desired  form 


A* 

g 


1 - (1/2)  tan  9 A9  -(1/6)(A9)2  ♦ (1/24)  tan  9 (A9)  3 ♦ ... 
tan  9 +(l/2)A0  - (1/6)  tan  0 (A0)2  -(l/24)(A0)3  + 
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To  cast  Eq.  2.5  into  a useful  form,  separate  it  into 


At  = i { In  [ -L— fP.l9-+--A0) 
g g 1 1 + sin  8 


] - In  [ 


cos  (9  + A9) 
cos  9 


] > 
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expand  the  arguments  in  power  series  in  A0,  and  then  make  use  of 


in  (1  + x)  = x - (1/2)  x2  + (l/3)  x3  -(l/4)x4  + ...  2.9 

which  is  valid  for  |x|  < 1.  Combination  of  these  steps  and  use  of  trigo- 
nometric identities  to  simplify  the  coefficients  in  the  resulting  power 
series  culminates  in 


At  - i [ 

g g L cos  0 


Ae  +^1/2)^FT  (A0)2  +^/6)- 


+ sin^  9 
cos3  0 


(A0) ' 


...].  2.10 


B.  Isospeed  Layer 

Refer  to  Fig.  1(b)  and  now  consider  the  properties  of  a ray  traversing 
a layer  in  which  the  speed  of  sound  cy  is  constant  over  its  thickness  Az. 

-^Recall  that  foT  both  models  the  same  value  of  the  ray  invariant  must  be 
used.  Thus,  we  assume  that  (cos  has  exactly  the  same  value  for  the 

straight-line  ray  in  the  isospeed  layer  as  for  the  circular  arc  in  the  iso- 
gradient layer.  As  we  have  formulated  the  problem  here,  this  is  equivalent 
to  relating  cy  to  both  the  speed  of  sound  c at  the  bottom  of  the  isogradient 
layer  and  the  gradient  g found  within  that  layer. > 

If  the  angle  of  elevation  is  0 + d0  for  the  ray  in  the  isospeed  layer, 
then  the  circular  arc  in  the  isogradient  layer  will  possess  the  same  slope 
at  that  depth  for  which  the  speed  of  sound  c(z)  is  equal  to  cy.  This  provides 
the  relation 

cy  ■ c(l  - | dz)  2.11 

where  dz  is  measured  from  the  bottom  of  the  layer.  Use  of  the  convenient 
quantity 
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h = & dz 
c 

and  application  of  Snell's  Law  yields 


cos  0 cos  (0  + d0) 
c = c (1  - h) 


from  which  we  extract  the  expression 

h = tan  0 d0  + (l/2)(d0)2  - (1/6) tan  0 (d0)3  -(l/24)(d0)4  + ...  2.14 


in  the  same  way  as  Eq.  2.3  was  obtained. 

As  an  aid  in  understanding  the  significance  of  h,  notice  that  if  h = 0, 
the  two  rays  have  the  same  slope  at  the  bottoms  of  the  layers,  and  if 
h = D/2,  the  rays  have  the  same  slope  at  the  middles  of  the  layers.  For 
this  latter  condition,  the  straight  ray  is  traveling  as  if  it  were  in  an 
isospeed  layer  possessing  speed  of  sound  c(l  - D/2),  which  is  the  average 
over  depth  of  the  speed  of  sound  in  the  isogradient  layer.  Notice  that, 
as  defined  above,  h must  always  be  less  than  or  equal  to  D. 

It  is  elementary  to  obtain  the  increment  in  range  A/tv  and  the  time-of- 
flight  Aty  for  the  ray  in  the  isospeed  layer: 


A*v  = cot  (0  + d0)  Az  = 


1 - tan  0 tan  d0 
tan  0 + tan  d0  2 


2.15 


and 


At  _ 1 Az 

v - c (1  - h)  sin  (9  + d0) 


2.16 
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Use  of  Eq.  2.13  and  the  definition  of  D allows  the  time-of-flight  to  be 
written  in  the  form 


1 _ & cos  (9  + d9) 

Aty  - D cos  6 


sin  (0  + d0) 


2.17 


Performing  the  same  kind  of  manipulations  as  before  (expansion  into  power 
series  in  the  incremental  angle,  use  of  Snell's  Law,  use  of  Eq.  2.3,  and 
so  forth)  allows  Eqs.  2.15  and  2.17  to  be  rewritten  in  the  forms 


LK  = 
v 


1 - tan  6 d9  -(1/3)  tan  6 (d8)3  + 
tan  0 + d0  + (1/3)  (d0)  3 + 


Az 


2.18 


and 


^ sin  9 [1  + (cot  9 - tan  9)  d9  - 2(d8)2  - (2/3)(cot  0 - tan  0)(d0)3  +... 
v 

2.19 


C.  Ratios  of  Incremental  Ranges  and  Times 

We  are  concerned  with  the  differences  between  incremental  ranges  and  times - 
of-flight  for  ray  paths  with  equal  ray  invariants  in  the  isogradient  and  iso- 
speed layers.  Investigation  of  these  small  differences  will  reveal  techniques 
to  minimize  the  discrepancies  between  the  range  and  time  results  of  t-i  two 
models . 

For  the  ratio  of  the  incremental  ranges,  direct  combination  of  Eqs.  2.7 
and  2.18  yields 
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MV  = 1 - tan8d8  - Cl/3)  tan8  (dO)  3 ♦ . . . 

A*g  tanB  + d9  + (l/3)  (d0)  3 ♦ ... 

x tan9  + (l/2)A9  -(2/3)  tan  e[fl/2)  A9j 2 - (l/3)[(l/2)  A@]  3 ♦ ...  ^2.20 

1 - tan0[(l/2)  A0]  -(2/$[(l/2)  A9]2  + (1/3)  tan0[(l/2)  A0]  3 + ... 

For  the  ratio  of  travel  times,  combine  Eqs.  2.10  and  2.19  and  then  use 

Eq.  2.3  (the  series  representation  for  D)  to  obtain 

Atg  = tan9A9  [1  + (l/2)  tan9A9  + (l/6)  (sec20  + tan29)(A9)2  + ...1 

Atv  tan0A0  + (l/2)  (A0)2  - (l/6)tan9  (A0) 3 -(l/24)(A0)4  + ...  2 21 

X[  1 + (cot0  - tan9)  d0  - 2(d9)2  > . ..]  . 
Study  of  these  two  ratios  reveals  that  they  depend  very  strongly  on  the 
angle  of  elevation  9 of  the  ray. 

As  an  aid  in  the  analysis  to  follow,  notice  that  for  all  angles  of 

elevation  Eq.  2.3  reveals  that  the  term  tan  9-A9  must  always  be  of  order 

D or  less,  and  Eq.  2.14  reveals  that  tan  9*d9  must  similarly  always  be  of 

order  h or  less.  On  the  other  hand,  the  terms  A0/(tan  0)  and  d0/(tan  0)  can 

be  arbitrarily  large  or  small,  depending  on  whether  0 is  small  or  large. 

As  a result,  it  will  often  be  necessary  to  investigate  these  ratios  for 

large  and  small  angles  separately.  The  distinction  between  large  and 

small  angles  of  elevation  must  be  made  on  the  basis  of  the  relationships 

between  the  angle  of  elevation  0 and  D or  h.  It  is  seen  for  Eqs.  2.3  and 

2.14  that  when  0 is  reduced  until  it  becomes  of  order  D (in  Eq.  2.3)  or 
}< 

of  order  h (in  Eq.  2.14)  then  the  incremental  angles  A0  and  d0  become  of 
the  same  respective  orders  and  the  terms  d0/£an  0)  and a0 /(tan  0)  change  from 
being  less  than  unity  to  being  greater  than  unity.  Since  typical  values 
of  D are  of  order  10"4,  it  is  clear  that  the  two  extremes  of  behavior 
occur  for  Jingles  0 much  less  than  or  much  greater  than  about  1°.  Thus, 
it  is  necessary  to  investigate  the  behaviors  of  Eqs.  2.20  and  2.21  for  the 
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cases  0 « h (since  h is  always  less  than  or  equal  to  D) , and  0 >>  D . 

In  what  follows,  we  must  distinguish  between  measuring  quantities  in 
terms  of  their  orders  in  A0  and  d0  and  measuring  quantities  in  terms  of 
their  order  in  D.  The  notation  we  adopt  to  make  this  distinction  is  as 
follows : 

(1)  If  a term  is  of  order  (A0)n  or  (d0)n,  it  will  be  represented  by 
0(n) . 

(2)  If  a term  is  or  order  Dn,  it  will  be  represented  by  0(Dn),  and 
similarly  for  h. 

The  expansions  in  incremental  angles  will  be  convenient  for  mathematical 
manipulations,  but  it  must  be  remembered  that  it  is  expansion  in  orders  of 
D which  will  reveal  the  sizes  of  the  discrepancies  between  ranges  and  times 
for  the  two  ray-tracing  models. 


D.  Solutionis  for  the  Incremental  Angles 

We  have  already  obtained  in  Eqs.  2.3  and  2.14  expressions  for  D and  h 
in  terms  of  the  incremental  angles  A0  and  d0  and  the  angle  of  elevation  0. 
Since  D and  h are  constants  for  the  layers,  these  expressions  are  implicit 
solutions  for  the  incremental  angles  in  terms  of  0 and  D or  h.  These 
expressions  must  be  inverted  to  obtain  the  desired  equations  for  A0  and 
d0. 

If  tan2  0 » D,  the  series  expansions  for  D and  h can  be  reverted  directly 
by  means  of  the  standard  formula: 

If 


a^x 


a2x 


a3x 


2.22 


1 1 -8 


then 


2 2.1 


X - -j  Y + C2  - aia3)  y" 


This  results  in  the  expressions 


D _ 1_  D 
tan  9 _ 2 . 3 


+ ... 


tan  9 


tan  9 


tan3  9 


2 2 

where  terms  just  through  0(D)  and  0(h)  have  been  retained. 

If  tan2$  <<  h , expressions  for  A 9 and  d9  can  be  obtained  in  the 
following  way:  Use  the  small  angle  approximation  tan  9=9  and  express 
the  incremental  angles  as  power  series  in  9 with  undetermined  coefficients , 
Substitute  these  into  Eqs.  2.3  and  2.14,  collect  terms  of  like  order,  and 
equate  to  zero.  This  reveals  that  the  assumed  power  series  are  valid  and 
evaluates  the  coefficients.  The  results  through  0(2)  are 


A9  = f~2 D - -6  + 9' 

2 /“2D1 


d9  = / 2h  ’ - 9 + — 9* 

2 /"inr 


These  expressions  are  useful  for  the  validation  of  some  equations  in  the  limit 
9-^0.  Otherwise,  they  are  of  no  practical  importance. 


1 1 -9 


If  A0  « 1 , then  0(3)  terms  can  be  dropped  in  Eqs.  2.3  and  2.14  and 
the  resulting  quadratic  equations  solved.  We  obtain 

A0  = / tan2  0 + 2D  - tan  0 2.28 

and 

d0  = / tan2  0 + 2h  - tan  0 . 2.29 

Note  that  for  the  situations  studied  here  (D  <<  1),  the  condition  A0  « 1 
will  be  met  for  all  0.  -<It  might  appear  that  these  approximate  solutions 
would  be  questionable  when  tan  0 gets  large.  However,  it  must  be  recalled 
that  tan  0‘A0  is  always  of  0(D)  or  less  which  means  that  the  discarded  0(3) 
term  will  always  be  0(D)  less  than  the  0(2)  term  in  each  of  Eqs.  2.3  and 
2.14.V  While  Eqs.  2.28  and  2.29  could  be  used  for  all  angles,  it  will  often 
be  more  convenient  to  use  the  simpler  forms  given  by  Eqs.  2.24  and  2.25  for 
large  angles  and  Eqs.  2.28  and  2.29  otherwise. 

This  last  case,  being  valid  for  all  0,  allows  evaluation  of  the  errors 
encountered  in  using  the  simpler  approximations  for  large  angles. 

If  we  assume  tan20  >>  D and  expand  Eq.  2.28  in  a power  series  in  D, 
it  is  seen  that  Eq.  2.28  and  Eq.  2.24  agree  identically  through  0(D2). 
Numerical  comparison  of  the  two  formulas  for  A0  reveals  that  they  agree  to 
within  0.5%  as  long  as  tan20  >_  10D,  with  the  discrepancy  decreasing  for 
larger  angles. 

If  just  the  first  term  is  retained  in  the  series  expansion  of  Eq.  2.24 
then  agreement  with  Eq.  2.28  is  within  0.5%  under  the  restriction 
tan20  >_  100D;  of  greater  significance  for  the  discussion  to  follow,  agreement 
is  within  5%  for  tan20  > 10D. 
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III.  RAY  PATH  MATCHING 


Two  rather  obvious  approaches  suggest  themselves  in  attempting  to 
reduce  the  differences  between  ray  paths  and  travel  times  in  isogradient 
and  isospeed  layers: 

(A)  Express  the  ratios  obtained  previously  in  terms  of  D and  h 
and  then  find  a relationship  between  D and  h which  tends  to  reduce  the 
differences  between  range  increments  and  between  times  of  flight. 

(B)  Expand  these  ratios  in  terms  of  the  incremental  angles  A0  and 
d0  and  find  a relationship  between  A0  and  d0  which  tends  to  reduce  the 
differences  between  range  increments  and  between  times  of  flight. 

It  might  appear  that  these  two  approaches  are  equivalent.  However, 

we  shall  see  that  there  will  be  distinct  differences  between  the  two 

when  the  angle  of  elevation  becomes  small. 

The  two  equations  to  be  investigated  are  Eq.  2.20 

1 - tan0  d0  - (l/3)tan0  (dB)  3 

Mg  1 - (l/2)tan0  A0  -(1/6)  (A0)  2 + (l/24)  tan0  (A9)  3 

x tan9  » (1/2)  A9  - (l/6)tanB (AB) 2 -(l/24)(A6)3 
tan0  + d9  + (l/3)(d9)3 

and  Eq.  2.21 

Atg  tan9  + (1  - tan20)  d9  - 2 tan9(d9)2 

Atv  tanB  + (l/2)A0  - (l/b^tan©  (A0) 2 2 

X[1  + (l/2)tan9  A0  +(1/6)  - Sin-  -9  (A0)2] 

COS20 

both  of  which  are  expressed  here  in  slightly  rearranged  forms. 

A.  Matching  Criterion  in  Terms  of  D and  h 

Because  the  expressions  relating  A0  and  d9  to  D and  h have  different 
forms  depending  on  the  values  of  the  angle  of  elevation,  it  is  necessary 
to  investigate  each  regime  separately. 


I j 
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(1)  Assume  tan29  >>  D 

For  this  case,  retain  terms  just  through  0(2)  in  Eq.  3.1,  revert  the 
series  in  the  denominator  with  the  help  of  Eqs.  2.22  and  2.23,  and 
collect  terms  of  like  powers  of  the  incremental  angles.  These  operations 
result  in 


AA 

v 

LK 

g 


1 . (l/2)A8  - de 
sin0cos0 


+ [4/4)  -(i/2). Ae  d9 ] . 

cos20  sin26  sin20cos2O  3.3 


Substitution  of  the  series  expansions  of  the  incremental  angles  in  terms 
of  D and  h (Eqs.  2.24  and  2.25)  and  collection  of  like-order  terms  in  D 
and  h gives 


AA 

g 


1 + (l/2)°  ~--2-  + (1/4)  [ 
sin20 


D2  + 6h2  - D2 

sin20  sin20tan20 


2Dh  j 
sin40 


3.4 


The  first  order  term  vanishes  for  h = D/2,  and  under  this  restriction  the 
ratio  of  incremental  ranges  becomes 


AA. 
v 

A fl 

g 


1 


I D2 

® tan20sin20 


tan20  > 10D . 


3.5 


The  condition  tan20  >_  10D  for  the  validity  of  Eq.  3.5  is  the  result  of 
retaining  just  through  0(2)  in  Eq.  3.3*.  When  (A0)2  and  (d0)2  are  calculated, 
terms  higher  than  0(2)  in  the  results  must  also  be  discarded,  but  this 
corresponds  to  retaining  just  the  first  terms  in  Eqs.  2.24  and  2.25.  Thus, 
the  term  of  0(D2)  in  Eq.  3.5  is  accurate  only  for  values  of  tan  0 such  that 
the  second  terms  in  Eqs.  2.24  and  2.25  are  negligible.  As  discussed  in 
the  previous  section  (after  Eq.  2.9),  this  requires  tan20  >_  10D  for  the  0(2) 
term  to  be  accurate  to  about  10%  of  its  value.  Analogous  manipulations  for 


l 


1 1 1 - 2 


the  ratio  of  incremental  times  (Eq.  3.2)  reveals  that  the  first  order 
term  also  vanishes  for  h = D/2  and  the  ratio  becomes 

At  1 r)2  i o 

— 8.  A 1 + ±-  — — ( — — + 4-  - 2 cos20),  tan26  > 10D.  3.6 

Atv  8 sin20  tan20  3 

(2)  All  angles,  and  tan20  « 1 

If  we  represent  the  angle  of  elevation  in  terms  of  D, 


tan20  = yD 


3.7 


where  y is  the  variable,  substitute  into  Eqs.  3.1  and  3.2,and  express  the 
incremental  angles  by  Eqs.  2.28  and  2.29,  the  ratios  of  incremental  ranges 
and  times  of  flight  become 

— » i — *-■  ~ [ 1 - yD (/  1 + 1/y'  - ✓ 1 + 2/y' /2  - 1/2)]  3.8 

/ 1 + 1/y' 

and 

Atg  = / 1 + 1/y'  - yDC3  - 3/  1 + 1/y'  + 2/y  1 

Atv  fl/2i(/  1 + 2/y'  + 1)  -(l/3)yD(l  - / 1 + 2/y’  + 1/y  ) 


X 


1 + (1/6)  yD(l  - / 1 + 2/y'  + 2/y) 

+ (2/3^ (yD) 2 (1  - / 1 + 2/y' ♦ 1/y) 


3.9 


1 1 1 - 3 


4 


If  we  assume  y » 1,  expansion  of  Eq.  3.8  yields 


vr  = 1 - J (l/y2  + D/y  - D2/3)  + 0(l/y3,  D/y2,  D2/y,  D3)  + .... 
g 

If  we  desire  about  10%  accuracy  in  the  second  term,  then  D/y  and  D2/3  can 
be  discarded  if  D/y  < 1/(10  y2) . This  gives  an  upper  limit  to  0 such 
that  tan20  _<  1/10.  Analysis  of  Eq.  3.9  is  similar  and  leads  to  about 
the  same  upper  limit  on  tan20.  Thus,  for  angles  satisfying  this 
inequality,  Eqs . 3.8  and  3.9  can  be  simplified  to 


A\ 

A*. 

g 


= (1/2) 


/ tan2@  + 2D ' + tan0 
/ tan20  + O' 


tan20  < 1/10 


3.10 


and 

At 

^ = (M- tan20  _<  1/10  . 3.11 

v g 

For  the  upper  bound  on  0 in  Eqs.  3.10  and  3.11  to  overlap  the  lower 
bound  on  0 in  Eqs.  3.5  and  3.6,  we  must  have  1/10  >_  10D  which  requires 
D _<  10'2. 

-(If  there  is  a range  of  angles  such  that  1/10  < tan20  <10D,  then  the 
full  forms  Eqs.  3.8  and  3.9  must  be  used  if  it  is  desired  to  estimate  the 
fractional  error  [1  - (A >tv/Mg)]  to  better  than  about  10%.  The  condition 
D ^ 10  2 would  be  virtually  impossible  to  encounter  when  accurate  ray 
tracing  is  being  done.> 

Some  pertinent  numerical  comparisons  are  given  in  Table  3.1  for  the 
ratios  of  incremental  ranges  predicted  by  Eqs.  3.5  and  3.10.  The  fractional 
errors  are  consistent  within  20%  for  tan20  >_  10D,  within  10%  for  tan20  > 20D, 


III -4 


h = D/2 


tan20 


Eq.  3.5 


Eq.  3.10 


1 C 

2 


- 1.25  x 10 


- 3.13  x 10 


- 5.00  x 10' 


- 1.25  x 10 


Comparison  of  ratios  of  incremental  ranges  for  the  case 
b = D/2  and  the  restriction  tan20  _<  1/10.  [The  ratios 
of  times-of-flight  (Eqs.  3.5  and  3.10)  deviate  from 
each  other  similarly.] 

Table  3.1 


1 1 1 -5 


and  within  1%  for  tan20  >_  100D.  Similar  calculations,  not  shown,  for  the 
ratios  of  incremental  times  reveal  comparable,  but  slightly  better,  agree- 
ment between  Eqs.  3.6  and  3.11  at  each  value  of  (itan2$/D. 

For  very  small  angles  the  fractional  deviations  are  quite  significant: 
From  Eqs.  3.10  and  3.11,  if  the  ray  is  horizontal  at  the  bottom  of  the 
isogradient  layer  then  the  ratio  of  incremental  ranges  approaches  1//  2 
and  the  ratio  of  travel  times  becomes  / 2.  Thus,  while  the  matching 
criterion  h = D/2  yields  very  satisfactory  agreement  for  the  incremental 
ranges  and  times  at  large  angles  of  elevation,  agreement  becomes  extremely 
poor  for  angles  approaching  grazing. 

B.  Matching  Criterion  in  Terms  of  A9  and  d8 

Direct  examination  of  Eq.  3.1  shows  that  if  we  set  d0  = A0/2,  all  terms 

of  0(1)  in  numerator  and  denominator  become  identical.  Now,  return  to  the 

fundamental  definitions  of  LA  and  LA  : Combination  of  Eqs.  2.4  and  2.15 

v g n 

with  d9  = A8/2  yields 


A* 

V _ COS0 1 

A-ig  “ tan (9  + A0/2)  sin(0  + A9)  - sin0 


gAz 

c 


Use  of  Eqs.  2.1  and  2.2  results  in 


cos(0  + A9/2)  cos8  - cos (9  + A8)  3.12 

LA  sin(0  + A9/2)  sin(0  + A9)  - sin9 
© 


Now, 


sin(9  + A9)  - sin0  = 2 sin(A9/2)cos(0  + A0/2) 
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and 


cos6  - cos (9  + A9)  = 2 sin(A9/2)sin(9  + A0/2) . 
Substitution  of  these  into  Eq.  3.12  shows  that 


^v 

A* 

g 


1 


3.13 


for  all  angles  of  elevation.  (It  is  worth  pointing  out  that  this  can  also 
be  seen  from  geometry:  In  the  layer  the  straight-line  ray  of  angle  of 
elevation  9 + A9/2  is  tangent  to  the  circular  arc,  subtending  angle  A9, 
at  its  midpoint,  and  is  therefore  parallel  to  the  chord  of  the  arc.) 

Direct  examination  of  Eq.  3.2  shows  that  the  criterion  d9  = A0/2 
similarly  causes  terms  through  0(1)  to  vanish.  Substitution  of  this 
choice  and  simplification  yields 

Atg  = tan9  -t-(l/2N)A9  + (1/12)  tan9  (tan29  - 1)(A9)2  3 H 

Atv  tan9  +(l/2)A0  - (l/6)tan9 (A9) 2 

Examination  of  the  limit  tan29  >>  D reveals 
At 

= 1 + (1/12)  (tan29  + 1)  (A0) 2 , tan20  >_  10D  3.15 

v 


where  the  restriction  tan29  >_  10D  follows  as  before.  Use  of  Eq.  2.24 
allows  the  above  expression  to  be  stated  in  terms  of  D, 


At 

At 


1 +('1/12) 


D2 

sin20 


tan29  > 10D  . 


3.16 


At  the  other  extreme,  under  the  restriction  tan29  <<  1,  Eq.  3.14  assumes 
the  approximate  form 
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tan20  < 1/10  . 


3.17 


Within  the  approximation  tan29  « 1,  for  which  tan0  = sin0,  Eqs . 3.18  and 
3.16  are  thus  equivalent.  Comparison  of  Eqs.  3.16  and  3.19  reveals  that 
the  fractional  error  between  travel  times  are  consistent  within  10%  if 
[1  - (tan9/sin0)2]  _<  1/10.  This  yields  the  upper  limit  on  tan20  stated  in 
Eqs.  3.17  - 3.19. 

Table  3.2  presents  some  values  of  Eqs.  3.16  and  3.18  as  functions  of 
(tan20)/D.  Notice  that  Atg/Aty  exhibits  a maximum.  Analysis  of  Eq.  3.18 
reveals  that  it  occurs  for  (tan20)/D  = 1/4  and  has  value  Atg/Atv  = 1 + D/24. 
For  D = 2 x 10’4  the  maximum  value  of  the  ratio  of  travel  times  is  thus 
1 + 8.33  x 10  6.  This  is  in  clear  contrast  with  the  results  of  the  case 
h = D/2  for  which  the  ratio  has  maximum  value  / 2 ' at  0 = 0. 
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t! 


de  = A9/2 


Eq.  3.16 


At 

At 

v 


Eq.  3.18 


tan29 

“D- 


1 + 


i2  sin29 


1 + 1_  tan9  (/  tan29  * r-taB-9I2 


/ tan29  + 2D  + tan9 


0 
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j 0.02 
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I 
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1 + 8.33  x 10"3  D 
1 + 4.17  x 10‘3  D 
1 + 1.67  x 10'3  D 


1 + 8.33  x 10  D 


1 + 0 


1 + 1.91  x 10  D 


1 + 2.47  x 10'2  D 
1 + 3.29  x 10'2  D 


1 + 3.83  x 10  D 


1 + 4.14  x 10"2  D 


1 + 3.93  x 10  D 


1 + 3.27  x 10"2  D 
1 + 7.25  x 10"3  D 
1 + 3.88  x 10"3  D 
1 + 1.62  x 10"3  D 


1 + 8.21  x 10  D 


Comparisons  of  ratios  of  incremental  times-of-f light 
for  the  case  d9  = A9/2  and  the  restriction  tan29  _<  1/10. 
The  ratio  of  incremental  ranges  is  exactly  unity  for  all 
angles.  v 

Table  3.2 
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IV.  TRACE -BACK  ERRORS 


Given  the  time-of-f light  of  the  received  signal,  the  value  of  the  ray 
invariant  over  the  path  it  has  travelled,  and  the  speed  of  sound  profile, 
the  received  signal  can  be  traced  back  to  its  apparent  source.  In  this 
section  we  will  treat  the  isogradient-layer  model  of  the  profile  as  our 
reference,  and  compare  the  discrepancies  in  range  and  depth  of  the  apparent 
location  of  the  source  between  isogradient-layer  models  and  isospeed-layer 
models  for  the  two  matching  criteria  h = D/2  and  d6  = A0/2. 

Let  us  define  the  quantity 


AA. 

1 


A A . 

A A . - AA  . = A*  ■ (1  - v/'1-  ) 
gi  vi  gi  ^ A A . ' 

gi 


4.1 


which  is  the  discrepancy  between  the  incremental  ranges  obtained  in  the  ith 
isospeed  and  isogradient  layers,  and  the  similar  quantity  for  the  difference 
between  the  times-of-flight  through  each  of  these  layers. 


At. 

l 


At 


vi 


At  .( 
vi 
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The  right  hand  side  of  each  of  these  equations  expresses  AA^  or  At^  in  terms 
of  the  ratios  obtained  in  Section  III. 

Examination  of  Figures  4.1(a)  and  4.1(b)  reveals  that  the  altitude  excess 
AZ  and  range  deficit  AR  for  the  source  position  from  trace -back  in  isospeed 
layers  (compared  to  that  for  trace-back  in  isogradient  layers)  are 

f 

AZ  = cc  sin  dr  .1,  At.  4.3 

f f i=l  l 
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.1 


Source  at 
position 

<W 


layer  model 


(a) 


AR  = - R£  = range  deficit 

(b) 


Figure  4.1 


and 


f 

AR  = Z AA.  - cot  9-  AZ 
i f 


4.4 


where  and  9^  are  the  speed  of  sound  and  angle  of  elevation  of  the  ray 
in  the  last  isospeed  layer. 

Since  we  have  usually  expressed  AA  ./An.  . and  At  ./At  . in  the  form 

r vi  gi  gi  vi 

[1  + (correction  term)^]  where  the  correction  term  is  found  to  reasonably 
high  accuracy,  we  can  in  what  follows  relax  the  rigor  maintained  up  to  this 
point,  since  we  are  concerned  more  with  the  approximate  sizes  of  AR  and  AZ 
rather  than  their  precise  values. 


A.  Trace-Back  Error  Formulas 
(1)  The  Case  h = D/2. 

(a)  Calculation  of  AA^: 

For  tan2  9 > 10D,  we  find  An.  ./An.  . from  Eq.  3.5,  and  A A . = Az.cot9. 

— vi  gi  n > gi  1 1 

from  Eqs.  3.5  and  2.15.  For  tan29^  _<  1/10,  the  ratio  comes  from  Eq.  3.10,  we 

can  replace  tan9.^  by  its  argument  9^^  with  negligible  error,  and  AAgi  comes 

from  combination  of  Eqs.  2.15  and  2.29,  AA  . = 2 Az./(/  9?  + 2D.*  + 9.). 

n gi  1^111 

Putting  these  together  gives 


r 


AA.  =» 
l 


i Di 
8 


AZ. 


tan39.sin29.  1 
l i 


if  tan29.  > 10D. 

l — x 


( — — ) Az.  if  tan29.  < 1/10 

/ 9?  + 2D.'  + 9.  / 9?  + D.'  1 1 


ill  li 
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(b)  Calculation  of  At.: 

1 

For  tan20.  > 10D. , the  ratio  is  found  from  Eq.  3.6  and  At  . = 
i—i  vi 

Azi/ (cisin®i^  results  from  Eq.  2.16.  For  tan20i  _<  1/10,  Eq.  3.11  yields 
the  ratio  and  since  angles  are  small  tan0^  can  be  replaced  with  0^; 

Eqs.  2.16  and  2.29  similarly  yield  Aty^  = Az^/ (c V 0?  + d!) . Thus  we  have 


r i Di 


o ( — + i-  - 2 cos20.)Az. 

8 sin30.  tan20 . 3 1 1 

l l 


( - )Az. 

v / 0?  + 2D.  + 0.  / 0*  + D.'  1 

ill  li 


if  tan20.  > 10D. 

l — l 


if  tan20.  < 1/10 

l — 
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where  simplification  has  been  achieved  by  the  additional  approximation 
c^  = c^  for  all  i.  (For  the  values  of  c typically  found  in  sea  water,  this 
approximation  is  good  to  within  a couple  per  cent.) 

(2)  The  Case  d0  = A0/2. 

We  have  seen  that  A/l^/AA.^  is  identically  unity  for  all  0^, 

A*.  = 0 all  0^  . 4.7 

Analysis  for  At^  proceeds  as  in  the  previous  case  except  that  Eqs. 
3.16  and  3.18  are  used  instead  of  Eqs.  3.6  and  3.11.  The  results  are 

r 

CfAtf  = J 

< 


1 

12 


D2 

1 — Az. 

sin30. 


if  tan20.  > 10D. 

l l 


T 9> 

3 l 


// 

\ 


0?  + 2D. 

l l 


0. 

l 


/ 0?  + 2D.'  + 0. 
ill 


Az. 

1 


if  tan20^  _<  1/10 
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The  formulas  for  these  cases  h = D/2  and  d0  = A9/2  are  summarized 


in  Figure  4.2. 


B.  Analysis  for  a Sample  Profile 

Assume  that  the  receiving  array  is  at  a depth  of  210  m,  the  source  is 
at  the  surface,  the  speed  of  sound  at  the  surface  is  1.5  x 10 3 m/sec,  Az 
is  3 m for  every  layer,  and  D = 2 x 10  4 for  all  layers.  This  corresponds 
to  a speed-of-sound  profile  which  is  nearly  isogradient  with  g = 0.1/sec. 

The  depth  deficits  Az  and  range  deficits  AR  of  the  source  location 
were  calculated  for  the  two  isospeed  models  as  functions  of  the  arrival 
angle  at  the  array.  The  calculations  used  the  approximate  equations  in 
Figure  4.2.  The  results  are  presented  in  Table  4.3  and  Figure  4.4. 

For  this  sample  profile,  the  condition  tan20  > 10  D is  equivalent  to 
0 ^2.6°  and  the  condition  tan20  1/10  is  equivalent  to  0 < 17*.  The 
pair  of  equations  for  A'u  and  the  pair  for  At^  thus  overlap  for  5®  <_  0 <_  10®, 
since  for  this  range  of  arrival  angles  all  angles  encountered  along  the  ray 
paths  lie  between  5®  and  14®.  In  this  range  of  overlap,  it  is  seen  that  the 
predicted  values  of  Az  and  AR  from  the  different  equations  for  each  of  the 
cases  h = D/2  and  d0  = A0/2  are  in  good  agreement.  The  poorest  agreement 
occurs  for  AR  in  the  case  h = D/2.  This  results  from  the  fact  that  for 
tan20^  > 10  D the  expression  for  AR  is  the  difference  of  two  large  nearly 
equal  quantities 

f f 

AR  * Z Ar^  - cos  0f  Z c _ At^ 

i=l  1 i=l  1 

each  of  which  is  correct  to  within  a few  percent.  The  disagreement  can  be 
seen  to  be  consistent  with  the  approximations  made  in  obtaining  the  large 
angle  equation  for  . As  0^  increases  the  discrepancies  in  AR  decrease 
rapidly  because  the  two  terms  become  more  dissimilar.  This  imprecision  in 
AR  does  not  occur  for  the  range  tan20^  < 1/10  because  in  this  case 
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D = 2 x 10  ’ 

Az  = 3 m 

Array  depth  = 210  m 
Source  at  the  surface 
c (surface)  = 1.5  x 103  m/sec 

Range  and  depth  deficits  for  the  two  isospeed  layer  models. 

Table  4.3 
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cf  At^  = Art^  so  that 

f 

Ar  = (1  - cos  9f)  E Art^ 
i=l 

and  the  difference  (1  - cos  0,)  does  not  involve  Art.  or  At. . 

f i i 

Figure  4.4  reveals  quite  clearly  that  the  use  of  the  isospeed  model 
based  on  h = D/2  results  in  very  large  errors  at  the  extreme  ranges. 

For  the  sample  profile  used  here,  it  is  seen  that  beyond  ranges  of  2000  m 
the  error  in  the  apparent  range  of  the  source  approaches  1 m and  the  error 
in  the  apparent  depth  exceeds  1 m. 

(It  is  to  be  noticed  that  the  errors  place  the  source  too  high  and  too 
near.  This  is  a result  only  of  the  model  used  for  traceback,  and  must  be 
distinguished  from  other  errors  which  may  enhance,  nullify,  or  reverse  this 
tendency . ) 
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V.  INITIALIZATION  ERROR 

A very  important  requirement  for  accurate  localization  of  the  source 
is  the  correct  value  for  the  ray  invariant.  An  error  in  this  quantity 
can  arise  from  incorrect  calculation  of  the  angle  of  arrival  0^  of  the 
sound  at  the  center  of  the  array  and  from  error  in  the  choise  of  the  value 
of  the  speed  of  sound  c^  at  the  center  of  the  array. 

Since  the  ray  invariant  is  cos  9r/cr  it  is  clear  that  the  effect  of 
an  incorrect  0^'  is  indistinguishable  from  that  of  an  incorrect  c\  Choosing 
.the  wrong  value  c^  is  equivalent  to  having  the  right  c^  but  the  wrong  0^  if 
c^  and  0^  are  related  by 

c 

cos  0'  = — r cos  0 5.1 

r c ' r 

r 

If  this  wrong  ray  is  traced  back  in  isogradient  water,  the  error  M.  becomes 


AK  = AA.  - AH' 

g g 


where  AV  is  for  the  incorrectly  chosen  ray.  Use  of  (2.4)  results  in 


M.  = 


g cos  0 


{ (sin  0 - sin  0^) 


cos  0 
r 

cos  0 ’ 
r 


(sin  0 ' - 


sin  0 ' ) } 
r 


5.2 


Analogously,  the  timing  error  is 

At  = At  - At' 

g g 

which,  with  (2.5)  yields 


At  * — Jin 

g 


1 + sin  0 
1 + sin  0 ' 


1 + sin  0 ' 

) 

1 + sin  0 

r 


5.3 
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If  we  assume  that  the  error  in  the  speed  of  sound  is  very  small,  so  that 


The  resultant  range  and  depth  deficits  are  then 


AZ  = (sin  9)  c At 
AR  = Art  - cot  0 AZ 

Since  both  At  and  Art  are  proportional  to  £,  we  see  that  the  errors  in 
locating  the  source  are  linear  in  the  fractional  error  £ in  the  speed 
of  sound  at  the  depth  of  the  center  of  the  array.  These  equations 
require  £ « tan2  0,  but  are  reasonably  accurate  as  long  as 
£/(tan2  0)  <_  1/4 . 

As  an  example,  let  us  assume  that  the  array  is  at  a depth  of  210  m, 
at  which  c^  = 1521  m/sec,  and  that  the  speed  of  sound  at  the  surface  is 
1500  m/sec.  This  isogradient  water  has  g = 0.1/sec  and  the  situation  is 
nearly  equivalent  to  that  analyzed  in  the  preceding  section.  For  a source 
at  the  surface,  the  range  and  depth  deficits  encountered  for  an  error  of 
0.15  m/s  in  the  speed  of  sound  at  the  array  (c^  = 1520.85  m/s)  are  pre- 
sented in  Figure  5.1. 

Notice  that  for  this  underestimate  of  c^  the  apparent  location  of  the 
source  is  significantly  too  deep  at  extreme  ranges.  For  the  sample  pro- 
file studied  here,  comparisons  of  these  errors  with  those  resulting  from 
the  isospeed  model  for  which  h = D/2  reveals  that  an  error  in  c^.  of  one 
part  in  104  is  far  more  serious  than  those  resulting  from  either  of  the 
isospeed  layer  models  for  the  speed  of  sound  profile. 
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9r 

R 

nominal 

AZ 

AR 

1° 

2270 

- 13  m 

- 2.2 

2 

2040 

- 5.8 

- 1.0 

5 

1520 

- 1.7 

- 0.32 

10 

1000 

- 0.56 

- 0.14 

20 

550 

- 0.15 

- 0.060 

50 

170 

- 0.014 

- 0.018 

c = 1521 
r 

c'  = 1520.85 
r 

g = 0.1 

Array  depth  = 210  m 
Source  at  the  surface 

Approximate  range  and  depth  errors. 
Table  5.1 
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VI.  RESULTS  AND  CONCLUSIONS 

The  analyses  of  the  preceding  sections  reveal  that: 

(1)  If  an  isogradient  layer  is  modeled  by  an  isospeed  layer  whose 
speed  of  sound  c^  is  the  average  over  depth  of  the  speed  of  sound  c(z) 
in  the  isogradient  layer. 


< c > 

g 


c tt  - f ) 
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the  trajectory  line  and  time  of  flight  through  the  layer  cam  be  seriously 
in  error  for  rays  whose  angles  of  elevation  are  small. 

(2)  If  in  the  isogradient  layer  the  ray  path  has  angle  of  elevation 
0 upon  entrance  to,  and  0 + A0  upon  exit  from,  the  layer,  then  if  the 
straight  line  ray  path  in  the  isospeed  layer  has  an  angle  of  elevation  of 
0 + (1/2)  A0,  the  error  in  the  trajectory  (at  entrance  and  exit)  vanishes 
identically  and  the  error  in  time  of  flight  remains  very  small  for  all 
angles  of  elevation. 

(3)  A very  small  error  in  either  0^  or  c^  can  introduce  serious  errors 
in  localization  of  the  source  when  the  ray  is  nearly  horizontal  over  some 
portion  of  its  path. 

Thus,  the  usual  method  of  basing  isospeed-layer  models  of  isogradient 
layers  on  the  average  over  depth  of  the  speed  of  sound  cy  = < c^  > can 
lead  to  serious  errors  in  the  localization  of  the  source  of  the  sound  if 
the  ray  becomes  nearly  horizontal  anywhere  on  the  path  between  source  and 
receiver.  If,  on  the  other  hand,  the  angle  of  elevation  of  the  ray  in  each 
isospeed  layer  is  required  to  be  0 + (1/2)  A0,  where  0 and  0 + A0  are  the 
entrance  and  exit  angles  of  the  ray  in  the  isogradient  layer,  then  the 
errors  become  acceptably  small.  It  is  worth  noting  that  this  approach  can 
be  shown  to  be  equivalent  to  requiring  c^  for  the  isospeed  layer  to  be 
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given  by 


c 

v 


< c > 

g 


, 1 (A0) 2 

8 i _ £ 

* c 

i -|  + ? (A0)2 

L 2 _ 

_ _ 
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It  is  our  opinion  that  better  than  either  of  these  choices  is  to  represent 
the  "true"  profile  by  isogradient  layers  for  which  the  speed  of  sound  is  allowed 
to  be  discontinuous  across  the  boundaries  of  the  layers.  This  method  of  approxi- 
mating the  profile  was  suggested  in  Figure  1.1(b).  It  is  to  be  noticed  that 
this  method  should  be  more  accurate  than  the  more  usually  encountered  con- 
strained fit  of  Figure  1.1(c). 

Finally,  very  small  errors  in  the  angle  of  elevation  or  speed  of  sound  at 
the  receiver  can  introduce  large  errors  in  localizing  the  source  if  the  ray 
is  nearly  horizontal  somewhere  over  its  path.  In  particular,  in  positive- 
gradient  water  the  most  important  part  of  the  speed-of-sound  profile  is  that 
closest  to  the  bottom  where  the  ray  is  most  nearly  horizontal . 
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